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Abstract 

A meta-surface made of a collection of nano-resonators characterized 
an electric dipole and a magnetic dipole was studied in the regime where 
the wavelength is large with respect to the size of the resonators. An 
effective description in terms of an impedance operator was derived. 


1 Introduction 

Metasurfaces are the bidimensional analogue of metamaterials [1]. They are 
made of resonant elements disposed on a surface. In this context, we study the 
field diffracted by a periodic set of linear nano-resonators, electromagnetically 
characterized by their scattering matrix S. We are interested in the regime 
where the wavelength is much larger than the size of the nano-resonators. Pos¬ 
sible with Mie-like resonances) or nano-wires doped with quantum dots. We 
proceed to an asymptotic analysis related to homogenization theory [1I1I1I5]. 
The field diffracted by the structure is derived and it is shown that it is charac¬ 
terized by an impedance operator. Our results extend to electric and magnetic 
dipoles the results in [B] where only electric dipoles were considered. 

2 Setting of the problem 

The structure under study is made of an infinite number of resonators invariant 
along z, periodically disposed at points Mp = {p x d,0), where d is the period 
and p G 7i. Each scatterer at position Mp is characterized in the frequency 
domain by a scattering matrix S{uj) as well as by electromagnetic parameters 
7 s and Ss- We assume that the wavelength in vacuum A = 2'kcIu: is much 
larger than the size of the resonators, which are assumed to be contained in 
a cylinder of diameter a. We therefore consider a linearly polarized incident 
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field where U'‘{x,y) = and a = fcosin0, 9 is an angle of 

incidence and ko = ^. According to the polarization 7 ^ and 6s can be either 
the relative permittivity or relative permeability. 

The stationary scattering problem considered therefore reads: Find a field 
U G satisfying div( 7 VU) + kg^U = 0 in the sense of Schwartz distri¬ 

bution in and such that = [/ — { 7 ® satisfies the outgoing wave conditions: 
y > 0,/TO(f7*VC/®) > 0 , 2 / < 0, Jm(C/®VC/®) < 0. The function 7 and 6 are 
equal, respectively, to 7 s and Ss inside the scatterers and to 1 outside. The 
following holds, see [7] 

Proposition 1. Apart possibly from a discrete set of wavenumbers ki, k 2 ,... 

The field U exists and is unique. 

Our point is to provide a simplified expression of the scattering problem by 
means of an impedance operator. 

Let Td denote the translation along x of amplitude d, ie. Td{f){x) = f{x—d). 

For later purpose, we define Y* = ]— K = a„ = a + nK and 
Pn = — For > fcg, we impose: i/3„ < 0. In the following, we denote 

T = R/dZ 

Let H = —div( 7 V-) — kgd and let us define for a &Y* the field of Hilbert 
spaces La(T) = G L^(T)} 

Lemma 1. The commutateur ofTd and H vanishes: [Td,H] = 0. 

Applying Floquet-Bloch analysis liiiiin], we obtain: 

Proposition 2. .■ The operator H has a direct integral decomposition H = 

fy, where Ha = —div( 7 V-)—kgd with domain D{Ha) = L^(T)) n H^^^{Y x R)} 


3 Multiple scattering approach 

The incident field has the expansion [TT]: U^{x,y) = E„a„J„(fcor)e™^ For one 
scatterer alone, the incident field gives rise to a field Up{x,y) = J2n ^nTn{x,y) 

where ipn{x,y) = Hn\kor)e'^'^^. For the infinite set of scatterers, this gives a 
diffracted field that reads as: 

U"{x,y) ='^sPipn{x -pd,y). (1) 

p,n 

Multiple scattering theory m allows to write that for p = 0: 

60 = (1 - SY)-^ Sa (2) 

where 6° = (..., 6° ^,..., 6°,.. .)^ and a = (..., a_„,..., a„,.. .)^. The 
matrix E is given by: 
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S(fco, ao) = Here proposition ([5]) was used through the intro- 

771^0 

ductionofaBlochphasee*“°'". In this expression, (Tom)^^ = Hp^2g{ko\m\d), 

that is 


Tr 


Om 


Ho{ko\m\d) -emHi{ko\m\d) i72(fco|w|d) 
emHi{ko\m\d) Ho{ko\m\d) -emHi{ko\m\d) 
H 2 {ko\m\d) emHi{ko\m\d) Ho{ko\m\d) 


where Cm = sign{m) (note that: = —sign(m). The following series |12) 

indexed by p appear: 


m^O 


and the entries of the matrix E(fco,Q;o) are: (S(A:o, Q;o))pg = ^p-q 

In the regime where k^a 1 the cylinder can be described by a 3 x 3 scat¬ 
tering matrix (this corresponds to an electric dipole and a magnetic dipole) and 
the field by 3 coefficients 6_i, b^, bi [U. Therefore, only 3 series are involved: 
Eq, Si, S 2 . It holds: 

/ Eo -Si S2 \ 

S(A:o,ao) = I Si Sq —Si 1 

\ ^2 Si So / 

In the extreme limit (a <C d) where the scatterers are very small as compared to 
the wavelength and the period, the scattering matrix S'(w) reduces to a scalar 
matrix so(w): the scatterers are thus dipoles with a dipole moment along ez and 
the only involved series is Sq, this situation was addressed in [^. The multiple 
scattering relation ([2]) then becomes: 

60(fco,ao) = (l-^oSo)“'^o. (4) 

where the series Sq can be written nniis]: 

So(A:o,ao) = ^ e*'='"^i7o(fco|m|d) (5) 

m^O 


For the more general case of an electric dipole and a magnetic dipole the fol¬ 
lowing asymptotic expressions hold in the limit k^d I m- 

Proposition 3. 


Eo(fco,ao)-1 - + fill 




(t) 


K 

TT^O 


( 6 ) 
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4 Scattering properties of the meta surface 


Define Pz = the electric moment and M = My) = ((6° + i{h\ — fo'Li)) 
the magnetic moment. We write m = + iMy and m* = — iMy and 

— {pm Pn)’! — {^nt Pn)- 

We can now state the following: 

Theorem 1. The total field has the expression: 

y> a: U{x, y) = e*(“o^-/5oy) + j.^^i(o.^x+is^y) 

y<a: C/{:r, y) = En 


(7) 


where: 

tn = d„o + (Pz + iM •«:„). 

Proof. We start with the following relation, obtained from Poisson formula: 


'rn = :^ {Pz + iM ■ n+), 

K 


( 8 ) 


^i/o(fco|r-nde.|)e*'="‘' = -^- 


,i(a„£c+/ 3 „|y|) 


Upon applying the operator d = dx + idy, using the fact that the series on the 
right hand side is normally convergent ( thanks to the term and using 

the relation: dHo{r) = —i/i(r)e®®, we obtain 


-koY,Ti{x-nd,y)P^-‘‘ = -J2 


2 {ictn Pm) iicr^X+Pr. 


|y|) 


(9) 


where e = sign{y). therefore we get: 

y) = 2 ^ b°o + ian{b°i + 6° i) + (6° 1 - b°)^ne ^,(^^x+0„ \v\)_ 


The result follows after some simple algebra. □ 

A simple, but interesting corollary is : 

Corollary 1. As n ^ +oo; 


K 2Kh\ K , 2A:6? 

r„ -- —m =-, ~ —m =- 

TT TT TT TT 


( 10 ) 


We are not a priori in the homogenization regime where kgd ^ 1 and hence 
there can be several reflected and transmitted orders. In expression 0, the 
propagative waves correspond to the fin’s that are real associated with the finite 
set U = {n € "Z, fin & R’*’} corresponding to the diffractive orders of the grating 
and the evanescent waves to the infinite set U'^ = {n € hfifin G K“}. 
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5 impedance operator formulation 


Our point is now to replace the set of nano-resonators by a meta-surface S which 
is simply the line y = 0. This requires to specify the boundary conditions there 
in terms of an impedance operator. 

To so, consider the continuation of the field U obtained by making a = 0 
in ©• The continued field is still denoted U. It is a singular distribution. To 
handle this situation, let us introduce the following fields of Sobolev spaces: 


UnC 


< +Oo| 

n J 


and the dual spaces: 

7J-/2(r)= L = + <+ool. 

n n J 

Let Z be the pseudo-differential operator defined by: Z[u\ = v where 
u{x) = ^ v{x) = ^ i^nUne"°‘"'" 


( 11 ) 


( 12 ) 


(13) 


It is straightforward to show the following: 

Proposition 4. Z is continuous and invertible from Ha *^^(T) to Ha ^0^); 
for s > 1. 

The inverse of Z is the admittance operator y defined by: z; = [m]. 


Let us denote F an element of Ha^^'^iY) x Ha^^^ ^(1^): representing the 
discontinuity of the field and its derivative through S. The traces of the field 

and F~ = — 

definition, it holds F+ -|- F~ = F. The Calderon projectors [T3] are defined by 
F+ = P+F, F~ = V~F. The preceding shows that 

Proposition 5. 


U{x,Q-) 

f(^>o-) 


• By 


and its derivative are: F^ = 


U{x,Q+) 

f(F0+) 


V+ = 


1 

2 


1 3^ 
Z 1 



(14) 


Obviously, it holds: + 'P~ = I, (^~)^ = V^V~ = 

= 0, as it should. 

The transmission conditions on § can be written: 


[U{x, 0+) - U{x, 0-)] = ^ ^ (15) 

TT ^ TT ^ iPnrn 


dU, du, 


= y] ^ , (16) 

„ „ t>«0 + rn 
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This suggests to define the following pseudo-differential operators, acting on 
u = 


T[u](a:) = ^ ( 

n 

1 + ^ \ 

^ OnO + *"n / 

(17) 

>V[u](a;) = 

2iKMy y 1 

(18) 


Both r„ and are bounded with respect to n (see eq. OT), hence the following: 
Proposition 6. The pseudo-differential operators X and W are isomorphism 

ofLim. 

These conditions can be rewritten conveniently in the operator form: 

Theorem 2. The traces of the field U{x,y) diffracted by the meta-surface under 
the illumination of a plane wave U’’ satisfy the impedance conditions: 

TF+ = F- , (19) 

where 

1 -yw 

-ZT 1 

The operator T is the transfer matrix of the meta surface. The discontinuity 
of the (effective) field 17 at y = 0 is due to the existence of a magnetic dipole 
moment. In the homogenization limit of very large wavelengths, i.e. larger than 
the wavelengths corresponding to magnetic resonances and larger than twice the 
period, there are only one transmitted and one reflected wave, the evanescent 
waves can be discarded and the magnetic resonances have no effect, we then 
have: 



Proposition 7. For a wavelength X > 2d and larger than the largest magnetic 
resonance, the propagative part of the field is given by: 


y>0: U{x,y) = + r^giio^o^+hov) 

y<0: U{x, y) = 

where ro = and t^ = l-\- r^. The transfer matrix becomes: 

r=[ ^ ^ 

' ‘ -WoTTr 1 


( 20 ) 

( 21 ) 


( 22 ) 


The form found for the transfer matrix matches that obtained in another 
context in [TO] . 
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6 Conclusion 


The scattering of a plane wave by a grating of nano-rods was described in the 
framework of meta-surfaces by exhibiting an impedance condition. This takes 
into account both the electric and the magnetic dipoles characterizing each 
nano-rod. This study can be generalized to higher multipoles but also to non¬ 
periodic, for instance quasi-periodic, structures to elementary scatterers 
deposited on an arbitrary smooth surface. A similar approach was used in m 
to study the coupling of a quantum emitter with the modes supported by the 
meta-surface, see also [16] . The proposed formalism can be used in this context 
to analyze the role of resonances. 
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